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Abstract 

Topological color codes are among the stabilizer codes with remarkable properties from quantum 
information perspective. In this paper we construct a lattice, the so called ruby lattice, with 
coordination number four governed by a 2-body Hamiltonian. In a particular regime of coupling 
constants, in a strong coupling limit, degenerate perturbation theory implies that the low energy 
spectrum of the model can be described by a many-body effective Hamiltonian, which encodes the 
color code as its ground state subspace. Ground state subspace corresponds to vortex-free sector. 
The gauge symmetry Z2 x Z2 of color code could already be realized by identifying three distinct 
plaqucttc operators on the ruby lattice. All plaqucttc operators commute with each other and with 
the Hamiltonian being integrals of motion. Plaquettes are extended to closed strings or string-net 
structures. Non-contractible closed strings winding the space commute with Hamiltonian but not 
always with each other. This gives rise to exact topological degeneracy of the model. Connection to 
2-colexes can be established via the coloring of the strings. We discuss it at the non-perturbative 
level. The particular structure of the 2-body Hamiltonian provides a fruitful interpretation in 
terms of mapping to bosons coupled to effective spins. We show that high energy excitations of 
the model have fermionic statistics. They form three families of high energy excitations each of 
one color. Furthermore, we show that they belong to a particular family of topological charges. 
The emergence of invisible charges related to the string-net structure of the model. The emerging 
fermions are coupled to nontrivial gauge fields. We show that for particular 2-colexes, the fermions 
can see the background fluxes in the ground state. Also, we use Jordan- Wigner transformation in 
order to test the integrability of the model via introducing of Majorana fermions. The four-valent 
structure of the lattice prevents the fermionized Hamiltonian to reduce to a quadratic form due 
to interacting gauge fields. We also propose another construction for 2-body Hamiltonian based 
on the connection between color codes and cluster states. The corresponding 2-body Hamiltonian 
encodes cluster state defined on a bipartite lattice as its low energy spectrum, and subsequent 
selective measurements give rise to the color code model. We discuss this latter approach along 
the construction based on the ruby lattice. 

PACS numbers: 03.65.Vf,75.10.Jm,71.10.Pm,05.30.Pr 
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I. INTRODUCTION 



Topological color codes (TCC) are a whole class of models that provide an instance of 
an interdisciplinary subject between Quantum Information and the physics of Quantum 
Many-Body Systems. 

Topological color codes were introduced class of topological quantum codes that 

allow a direct implementation of the Clifford group of quantum gates suitable for entangle- 
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ment distillation, teleportation and fault-tolerant quantum computation. They are defined 
on certain types of 2D spatial lattices. They were extended to 3D lattices [2| in order to 
achieve universal quantum computation with TCCs. This proposal of topological quantum 
computation relies solely on the topological properties of the ground state sector of certain 
lattice Hamiltonians, without resorting to braiding of quasiparticle excitations. In addition 
to these applications in Quantum Information, topological color codes have also a natural 
application in strongly correlated systems of condensed matter with topological orders. In 
3| was found that TCCs can be extended to arbitrary dimensions, giving rise to topological 
orders in any dimension, not just 2D. This is accomplished through the notion of D-colexes, 
which are a class of lattices with certain properties where quantum lattice Hamiltonians are 
defined. This corresponds to a new class of exact models in D=3 and higher dimensions 
that exhibit new mechanisms for topological order: i/ brane-net condensation; ii/ existence 
of branyons; iii/ higher ground-state degeneracy than other codes; iv/ different topological 
phases for D > 4 etc. In what follows, we shall focus only on 2D lattices. 

Physically, TCCs are exotic quantum states of matter with novel properties. They are 
useful for implementing topological quantum computation, but they have also an intrinsic 
interest by their own. Then, a natural question arises as to how to implement experimen- 
tally these new quantum systems by means of light, atoms or some other means. This is a 
challenge since TCCs are formulated in terms of Hamiltonians with many-body terms, the 
simplest having 6-body interactions in a hexagonal lattice. But the most common interac- 
tions in nature are typically 2-body interactions. 

There are several approaches to trying to solve this challenge, depending on the type of 
scenario we envision to be in practice and the practical rules we are supposed to be allowed 
to have at our disposal. 

Let us start first for what we may call a 'quantum control scenario'. By this we simply 
mean that we are able to perform very controllable quantum operations on our system 
that we have prepared artificially. In particular, we suppose that we can perform quantum 
measurements on the qubits and having ancilla qubits at will. Under these circums ;ances, we 
can resort to cluster states [4] and measurement-based quantum computation js, O]. This 
is because TCCs can be described by a certain cluster state construction [7] within this 
scenario. Then, it is possible to use a technique to obtain graph states as ground states of 
2-body qubit Hamiltonians js, 9|. We show this construction in |Al However, this scenario 
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is experimentally very demanding and it is left for the future when will it be achieved 
completely. Therefore, it is convenient to seek other alternatives. 

Thus, let us move onto a 'condensed matter scenario'. The terminology is intended just 
to be illustrative, rather than precise. In fact, the scenario goes beyond condensed matter 
and may well be a quantum simulation of our system by means of engineering a set of 
photons, atoms or the like. The important difference now is that external measurements on 
the system, or ancilla qubits, are not allowed in order to obtain the desired Hamiltonian for 
the TCCs. We want to remain in a framework based on Hamiltonians with solely 2-body 



interactions 10|. 



We have introduced a new quantum 2-body Hamiltonian on a 2D lattice with results 
that follow the twofold motivation concerning topics in Quantum Information and Quantum 
Many-Body Systems: 



i/ to achieve scalable quantum computation llMlSj: 



ii/ to perform quantum simulations with light, atoms and similar available means [l4j-|23|. 
This is so because, on one hand, the Hamiltonian system that we introduce is able to 

n 

reproduce the quantum computational properties of the topological color codes (TCC) pj- 
at a non-pertubative level as explained in Sect JI VI This is an important step towards 
obtaining topological protection against decoherence in the quest for scalability. On the 
other hand, the fact that the interactions in the Hamiltonian appear as 2-body spin (or 
qubit) terms makes it more suitable for its realization by means of a quantum simulation 
based on available physical proposal with light and atoms. 

In a framework of strongly correlated systems in Quantum Many-Body Systems, one of 
the several reasons for being interested in the experimental implementation of this Hamil- 
tonian system is because it exhibits exotic quantum phases of matter known as topological 



orders, some of its distinctive features being the existence of anyons |24j-|26|. In our everyday 
3D world, we only deal with fermions and bosons. Thus, exchanging twice a pair of particles 
is a topologically trivial operation. In 2D this is no longer true, and particles with other 
statistics are possible: anyons. When the difference is just a phase, the anyons are called 



abelian. Anyons are a signature of topological order (TO) |27|, l28|, and there are others as 
well: 

• there is an energy gap between the ground state and the excitations; 
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• topological degeneracy of the ground state subspace (GS); 

• this degeneracy cannot be lifted by local perturbations; 

• localized quasiparticles as excited states: anyons; 

• edge states; 

etc. 

These features reflect the topological nature of the system. In addition, a signature of 
the TO is the dependence of that degeneracy on topological invariants of the lattice where 
the system is defined, like Betti numbers 

But where do we find topological orders? These quantum phases of matter are difficult 
to find. If we are lucky, we may find them on existing physical systems such as the quantum 
Hall effect. But we can also engineer suitable quantum Hamiltonian models, e.g., using polar 



molecules on optical lattices IJ, ll5|, |20|, or by some other means. There are methods for 



demonstrating topological order without resorting to interfero metric techniques [29 1. 

In this paper we present new results concerning the realization of 2-body Hamiltonians 
using cluster states techniques on one hand, and without measurement-based computations 
on the other. In this latter case, we present a detailed study of the set of integrals of motion 
(lOM) in a 2-body Hamiltoinan, fermionic mappings of the original spin Hamiltonian that 
give information about the physics of the system and which complements previous results 
using bosonic mapping techniques jlO |. 

This paper is organized as follows: in Sect. [Til we present color codes as instances of 
topological stabilizer codes with Hamiltonians based on many-body interacting terms and 
then introduce the quantum Hamiltonian model based solely on 2-body interactions between 
spin-i particles. The lattice is two-dimensional and has coordination number 4, instead of 
the usual 3 for the Kitaev model. It is shown in Fig. |3]and it is called ruby lattice. In Sect. 
mil we describe the structure of the set of exact integrals of motion (lOM) of the 2-Body 
model. We give a set of diagrammatic local rules that are the building blocks to construct 
arbitrary lOMs. These include colored strings and string-nets constant of motion, which is a 
distinctive feature with respect to the Kitaev's model. In Sect. llVt we establish a connection 
between the original topological color code and the new 2-Body color model. This is done 
firstly at a non-perturbative level using the colored string integrals of motion that are related 
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with the corresponding strings in the TCC. Then, using degenerate perturbation theory in 
the Green function formahsm, it is possible to describe a gapped phase of the 2-Body color 
model that corresponds precisely to the topological color code. In Sect. IVj we introduce 
a mapping from the original spin-| degrees of freedom onto bosonic degrees of freedom in 
the form of hard-core bosons which also carry a pseudospin. This provides an alternative 
way to perform perturbation theory and obtain the gapped phase corresponding to the 
TCC. It also provides a nice description of low energy properties of the 2-Body model and 
its quasiparticles. In Sect. IVIt we introduce another mapping based on spinless fermions 
which is helpful to understand the structure of the 2-Body Hamiltonian and the presence 
of interacting terms which are related to the existence of stringnets constants of motion. 
Sect. IVIIl is devoted to conclusions and future prospects. |X] describes how to obtain 2-Body 
Hamiltonians for topological color codes based on cluster states and measurements using 
ancilla qubits. 



II. QUANTUM LATTICE HAMILTONIAN WITH TWO-BODY INTERACTIONS 



A. Color Codes as Topological Stabilizers 



Some of the simplest quantum Hamiltonian models with topological order can be obtained 
rom a formalism based on the local stabilizer codes borrowed from quantum error correction 



30| in quantum information 3l|, l32|. These are spin-| local models of the form 

H = -J2S., ^.eP„:=(^,a?^, <...,<,<). (1) 

i 

where the stabilizer operators 5*^ constitute an abelian subgroup of the Pauli group P„ of n 
qubits, generated by the Pauli matrices not containing —1. The ground state is a stabilizer 
code since it satisfies the condition 



Si\GS) = |GS),Vz, 



(2) 



and the excited states of H are gapped, and correspond to error syndromes from the quantum 
information perspective 

'5.1^) = -!^). (3) 
The seminal example of topological stabilizer codes is the toric code 



basically two types of known topological stabilizer codes [l3|. It is possible to study this type 



331]. There are 



of homological error correcting codes in many different situations and perform comparative 



studies 
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-42!. Topological color code (TCC) is the another relevant example of topological 
stabilizer codes, with enhanced computational capabilities jl-3|. In particular, they allow 
the transversal implementation of Clifford quantum operations. The simplest lattice to 
construct them is a honeycomb lattice A shown in FigJH where we place a spin-i system at 
each vertex. There are two stabilizer operators per plaquette: 



'1 '2 '3 '4 '5 



(4) 
(5) 



/ 

where r^'s = x,y) are usual Pauli operators. There exist six kinds of basic excitations. 




FIG. 1: The hexagonal lattice is an example of 3-colorable lattice by faces, and also by edges. A 
topological color code can be defined on it by associating two stabilizer operators for each plaquette 

To label them, we first label the plaquettes with three colors: Notice that the lattice is 
3-va,e,. a.d has 3-co,o.able p.a,„ettes. We call .,ch lattices 2-eo.exes g One ea. define 

color codes in any 2-colex embedded in an arbitrary surface. There exists a total of 15 
nontrivial topological charges as follows. The excitation at a plaquette arises because of the 
violation of the stabilizer condition as in (j3]). Consider a rotation applied to a certain 
qubit. Since anticommutes with plaquette operators Bj of neighboring plaquettes, it will 
put an excitation at corresponding plaquette. Similarly, if we perform a rotation on a 
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qubit, the plaquette operators are violated. These are the basic excitations, two types 
of excitations per each colored plaquette. Regarding the color and type of basic excitations, 
different emerging excitations can be combined. The whole spectrum of excitations is shown 
in FiglJI^a). Every single excitation is boson by itself as well as the combination of two basic 
excitations with the same color. They form nine bosons. However, excitations of different 
color and type have semionic mutual statistics as in Figl2t^b). The excitations of different 
color and type can also be combined. They form two families of fermions. Each family 
of fermions is closed under fusion, and fermions from different families have trivial mutual 
statistics. This latter property is very promising and will be the source of invisible charges 
as we will discuss in SectJVl The anyonic charge sectors are in one to one correspondence 
with the irreducible representations (irreps) of the underlying gauge group, and the fusion 
corresponds to decomposition of the tensor product of irreps. 



Bosons: 



Fermions 1: 
Fermions 2: 



y xy 



(a) 

/(•.•)=/{*,^)=/(^,»)=-1 

/(•.•)=/(^,«)=/(»,3')=-1 

(b) 



FIG. 2: (a) Classification of excitations for the topological color code model ([5]), nine bosons and 
two families of fermions (b) The nontrivial phase arising from the braiding of different charges. 



We describe all above excitations in terms of representation of the gauge group of the 
TCC. Before that, let us make a convention for colors which will be useful for subsequent 
discussions. We refer to colors by a bar operation c that transform colors cyclically as r = g. 
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g = b and b = r. The elements of the gauge group Z2 x Z2 are {e, r,b,g}. Each excitation 
carries a topological charge. The corresponding topological charge can be labeled by the pair 
(g, x), where g G Z2 x Z2 and x an irrep of this group [10|]. We label them as Xe(c) = Xc(c) = 
— Xc(c) = 1. Therefore, there are nine bosons labeled by (c, Xe), (e, Xc) and (c, Xc) and six 
fermions (c, Xc) and (c, xz)- Taking into account the vacuum with trivial charge (e, Xe), color 
code has sixteen topological charges or superselection sectors. Regarding the fusion process, 
fusion of two charges (g, Xc) and (g',Xc') give rises to {qq',XcXc') charge. Additionally, 
the braiding of charge (g, Xc) around charge {q',Xc') produces the phase Xc{q')Xc'{q)- An 
excitation at a c-plaquette has (c, Xe) charge if —B^ = = 1, {e, Xc) charge if B^ = —B^ = 
1 and (c, Xc) charge if B^ = B^ = —1. 

It is also possible to use both types of topological stabilizer codes, either toric codes or 
color codes, to go beyond homological operations. This corresponds to performing certain 
types of operations called code deformations, which may alter the topology of the surface 
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47|. 



allowing an extension of the computational capabilities of these 2D codes 

Active error correction procedures are particularly interesting in the case of topological 
stabilizer codes. They give rise to connections with random statistical mechanical models 
like the random bond Ising model for the toric code |47] and new random 3-body Ising models 
for color codes {3]. The whole phase diagram p — T has been mapped out using Monte 
Carlo, which in particular gives the value of the error threshold Pc- This particular point 
can also be addressed using multicritical methods 49[ . There is experimental realization 
of topological error correction [50]. Without external active error correction, the effect of 
thermal noise is the most challenging problem in toric codes 
effects of topological order in color codes has also been studied js? ]. 

In all, the type of entanglement exhibited by topological color codes is very remarkable 



58N64|. A very illustrative way to see this is using the connection of the ground state of 



topological codes with standard statistical models by means of projective measurements 
0, IgsItI]. For TCCs, this mapping yields the partition function of a 3-body classical 
Ising model on triangular lattices 7|. This 3-body model is the same found in active er- 
ror correcting techniques 48|], but without randomness since there is no noise produced by 
external errors. This type of statistical mapping allows us to test that different computa- 
tional capabilities of color codes correspond to qualitatively different universality classes of 
their associated classical spin models. Furthermore, generalizing these statistical mechan- 
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ical models for arbitrary inhomogeneous and complex couplings, it is possible to study a 
measurement-based quantum computation with a color code state and we find that their 
classical simulatability remains an open problem. This is in sharp contrast with toric codes 
which are classically simulable within this type of scheme 



B. The Model 

In nature, we find that interactions are usually 2-body interactions. This is because 
interactions between particles are mediated by exchange bosons that carry the interactions 
(electromagnetic, phononic, etc.) between two particles. 

The problem that arises is that for topological models, like the toric codes and color codes, 
their Hamiltonians have many-body terms ([5]). This could only achieved by finding some 
exotic quantum phase of nature, like FQHE, or by artificially enginering them somehow. 

Here, we shall follow another route: try to find a 2-body Hamiltonian on a certain 2D 
lattice such that it exhibits the type of topological order found in toric codes and color 
codes. In this way, their physical implementation looks more accessible. 

In fact, Kitaev 



72| introduced a 2-body model in the honeycomb lattice that gives rise 
to an effective toric code model in one of its phases. It is a 2-body spin-| model in a 
honeycomb lattice with one spin per vertex, and simulations based on optical lattices have 



been proposed 



17|. 



The model features plaquette and strings constants of motion. Furthermore, it is exa ctly 



solvable, a property that is related to the 3- valency of the lattice where it is defined j72|-l78|. 
It shows emerging free fermions in the honeycomb lattice. If a magnetic field is added, 
it contains a non-abelian topological phase (although not enough for universal quantum 
computation). Interestingly enough, another regime of the model gives rise to a 4-body 
model, which is precisely an effective toric code model. A natural question arises: Can we 
get something similar for color codes? We give a positive answer in what follows. 

Motivated by these physical considerations related to a typical scenario in quantum many- 
body physics, either condensed matter, AMO physics or the like, we will seek a quantum 
spin Hamiltonian with the following properties: 
i/ One of its phases must be the TCC. 

ii/ To have two sets of plaquette operators generating a Z2 x Z2 local, i.e. gauge, symmetry. 
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iii/ To have string-nets and colored strings lOM as in the TCC, but in all coupling regimes. 

Thus, the reasons behind demanding these properties are to guarantee that the sought 
after model will host the TCC. For instance, property i/ means that we must be able to 
generate the 2D color code Hamiltonian consistently at some lowest order in perturbation 
theory (PT). This we shall see in Sect JIVBl Likewise, properties ii/ and iii/ are demanded 
in order to have the fundamental signatures regarding gauge symmetry and constants of 
motions associated with TCCs. Notice that we have not demanded that the model be 
exactly solvable. This is a mathematical requisite, rather than physical. We leave the door 
open for considering larger classes of models beyond exactly solvable models, which may 
be very interesting and contain new physics. For example, according to those properties, it 
would be possible to have models with a number of lOMs that scales linearly with N, the 
number of spins or qubits. Thus, the Kitaev model has a number of lOMs of |iV. 

Our purpose is to present first the 2-body quantum Hamiltonian in 2D [io|, and then to 
analyze diverse possible mappings in later sections, like using bosonic and fermionic degrees 
of freedom. The analysis of the set of lOMs will play also a crucial role in the understanding 
of our model as we shall see in Sect illll 

It is a 2-body spin-1/2 model in a 'ruby' lattice as shown in FigJS] We place one spin per 
vertex. Links come in 3 colors, each color representing a different interaction. 



X, red links 



y, green links (6) 
z, blue links 

For a suitable coupling regime, this model gives rise to an effective color code model. 
Furthermore, it exhibits new features, many of them not present in honeycomb- like models: 

• Exact topological degeneracy in all coupling regimes (4^ for genus g surfaces). 

• String-net integrals of motion. 

• Emergence of 3 families of strongly interacting fermions with semionic mutual statis- 
tics. 

• Z2 X Z2 gauge symmetry. Each family of fermions sees a different Z2 gauge subgroup. 



12 



FIG. 3: A lattice with coordination number 4 where the 2-body quantum lattice Hamiltonian for 
the color codes is defined according to spin-spin interactions coded by the colors of the links, as in 
([6]). A plaquette can be distinguished by an inner hexagon, an outer hexagon and six blue triangles 
between them. 

III. STRING OPERATORS AND INTEGRALS OF MOTION 

We can construct integrals of motion (lOM), / G P„, [Hcc, I] = 0, following a pattern of 
rules assigned to the vertices of the lattice, as shown in FigJU These rules are constructed 
to attach a Pauli operator of type erf, af or erf to each of the vertices i. The hues around the 
vertices, either wavy lines or direct lines, are pictured in order to join them along paths of 
vertices in the lattice that will ultimately translate into products of Pauli operators, which 
will become lOMs. Clearly, af operators are distinguished from the rest. The contribution 
of each qubit in the string operator is determined in terms of how it appears in the string. 
Its contribution may be determined by the outgoing red and green links which have the 
qubit as their end point in the string. In this case the or cr^ Pauli operators contribute 
in the string lOM . If a typical qubit crossed only by a wavy line as shown in FigJU^a), it 
contributes a Pauli operator in the string. To have a clear picture of string operators, a 
typical example has also been shown in Fig|l](b). Part of string is shown on the left and its 
expression will be the product of Pauli operators which have been inserted in open circles 
on the right. With such definitions for string operators and their supports on the lattice, 
now we turn on to analyze the relevance of strings to the model. In particular, we will 
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construct elementary string operators with the local symmetry of the model. Therefore, in 
this way we are representing the local structure of the lOMs of our 2-body Hamiltonian 
([5]). We will illustrate them with several examples of increasing complexity. The ground 



FIG. 4: (a) A diagrammatic representation of the local structure of the integrals of motion of 
the 2-body Hamiltonian ([6]). The colored links represent different spin-spin interactions, (b) An 
example of contribution of Pauli operators in a string. 

state of a lattice model described by the Hamiltonian ([S]) is a superposition of all closed 
colored strings. Indeed, it is invariant under any deformations of colored strings as well as 
splitting of a colored string into other colors. In other words, the ground state is a string-net 
condensed state and supports topological order. The gauge group related to this topological 
order is Z2 x Z2. Such symmetry of topological color code can be realized via defining a 
set of closed string operators on the ruby lattice. We shall verify the gauge symmetry by 
identifying a set of string operators on the lattice of FigJHl 

Let us start by constructing the elementary string lOM as shown in FigjSl They are 
denoted as I = A, B, C. They are closed since they have not endpoints left. The elementary 
closed strings are plaquettes. By a plaquette we mean an inner hexagon and an outer hexagon 
with six triangles in between. For a given plaquette it is possible to attach three string 
operators. For each closed string, the contribution of Pauli operators are determined based 
on outgoing red and green links or wavy lines as in FigJH Let Vf stand for a set of qubits 
on a plaquette. Note that each plaquette contains 18 qubits corresponding to six triangles 
around it. For first plaquette operator in Figj5]we can write its explicit expression in terms 





(a) 



(b) 
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of Pauli matrices as 5*^^ = IliGy/ where / denotes the plaquette and u = x,y, depending 
on outgoing red or green hnks, respectively. Similarly the second plaquette operator has an 
expression as Sj^ = Yl^^y^ cTi- The third string is only a closed wavy string which coincides to 
the inner hexagon of the plaquette. It's expression is S'j^ = YlieVh ^t^ where Vh stands for six 
qubits on the inner hexagon. The three closed strings described above are not independent. 
Using the Pauli algebra, it is immediate to check that they satisfy S"^ = —SfSj^. Thus, 
there exist 2 independent lOMs per plaquette: this is the Z2 x Z2 local symmetry of the 
model Hamiltonian ([6]). 




A B C^-AB 

FIG. 5: Schematic drawing of the plaquette lOMs according to the local rules in FigHl There are 3 
lOMs denoted as A,B,C, but only 2 of them are independent. This corresponds to the symmetry 
Z2 X Z2 of the model. 

Plaquette operators commute with each other and with any other lOM. If a lOM corre- 
sponds to a nontrivial cycle c, it is possible to find another lOM that anticommutes with 
it, namely one that corresponds to a cycle that crosses once c' . Thus, lOMs obtained from 
nontrivial cycles are not products of plaquette operators. 

Each string operator squares identity since we are working with qubits. Plaquette oper- 
ators corresponding to different plaquettes commute with each other and also with terms in 
Hamiltonian in ([6]) since they share in zero or even number of qubits. Therefore, the closed 
strings with the underlying symmetry obtained above define a set of integrals of motion. 
The number of integrals of motions is exponentially increasing. Let 3A^ be the total number 
of qubits, so the number of plaquettes will be Regarding to the gauge symmetry of 
the model, the number of independent plaquette operators is A^. This implies that there 
are 2^ integrals of motion and allow us to divide the Hilbert space into 2^ sectors being 
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eigenspaces of plaquette operators. However, for closed manifold, for example a torus, all 
plaquette operators can not be independently set to +1 or —1 because they are subject to 
some constraint. All other closed string operators that are homologous to zero, i.e. they are 
homotopic to the boundary of a plaquette, are just the product of these elementary plaque- 
tte operators. It is natural that all of them are topologically equivalent up to a deformation 
and commute with the Hamiltonian of the model. 




FIG. 6: An example of a stringnet lOM. Notice the presence of branching points located around 
blue triangles of the lattice. This is a remarkable difference with respect to honeycomb models like 
the Kitaev model. 



The most general configuration that we may have is shown in Figj6l We call them string- 
nets lOM since in the context of our model, ^ they can be thought of as the string-nets 
introduced to characterize topological orders [79]. The key feature of these lOMs is the 
presence of branching points located at the blue triangles of the lattice. This is remarkable 
and it is absent in honeycomb 2-body models like the Kitaev model. When the string- 
nets lOM are defined on a simply connected piece of lattice they are products of plaquette 
operators. More generally, they can be topologically non-trivial and independent of plaquette 
operators. 

As a special case of lOMs we have string configurations, i.e., paths without branching 
points. They correspond to the different homology classes of the manifold where the lattice 
is embedded , and are needed for characterization of the ground state manifold. Some 
examples are shown in Figd They may be open or closed, depending on whether they 
have endpoints or not, respectively. Strings lOM are easier to analyze. String-nets lOM are 
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products of strings lOM. For a given path, there exist 3 different strings lOM. These are 
denoted as A, B, C in FiglTl We must introduce generators for the homology classes defining 




FIG. 7: Examples of standard string configurations of lOMs, i.e., without branching points. For 
each path, we can in principle make 3 different assignments of lOMs, but again only 2 of them are 
independent as with plaquette lOMs. This is another manifestation of the Z2 x Z2 symmetry of 
the model. 

closed manifold. Homology classes of the torus are determined by realizing two nontrivial 
loops winding around the torus. In the Kitaev's model there are only two independent such 
nontrivial closed loops. However, the specific construction of the lattice and contribution of 
the color make it possible to define for each homology class of the torus two independent 
nontrivial loops. These closed strings are no longer combination of plaquettes defined above. 
Let stand for such string, where A and /i denote the type and homology class of the 
string. For each homology class of the manifold we can realize three different types of string 
operators depending on how the vertices of the lattice are crossed by the underlying string. 
Each qubit crossed by the string contributes a Pauli operator according to the rules in FigJH 
Again, using Pauli algebra we can see that only two of them are independent, as with the 
plaquette lOMs. 
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where t is the number of triangles on the string. To distinguish properly the three types we 
have to color the lattice. We could already use the colors to label strings. Strings are then 
red, green or blue. This is closely related to the topological color code Q. The latter 
relation shows that each string can be constructed of two other homologous ones, which is 
exactly the expression of the Z2 x Z2 gauge symmetry. Each non-contractible closed string 
operator of any homology commutes with all plaquette operators and with terms appearing 
in the Hamiltonian, so they are constant of motions. But, they don't always commute with 
each other. In fact, if the strings cross once then 

K s:] = 0, (8) 

but 

' s;, St \ = 0. (9) 



This latter anticommutation relation is a source of exact topological degeneracy 80|] of the 
model independent of phase we are analyzing it. 



IV. A GAPPED PHASE: THE TOPOLOGICAL COLOR CODE 



A. Non-Perturbative Picture 



In this subsection we discuss the ruby lattice is connected to the 2-colex even at the 
non-perturbative level. Then, in the subsequent sections we verify it using quantitative 
methods. From the previous discussion on lOMs, we have already seen a connection with 
the topological color codes. Now, we want to see how different strings introduced above are 
related to coloring of the lattice. To this end, consider the closed strings A,B,C in FigJS] 
The closed strings A and B can be visualized as a set of red and green links, respectively. 
With such visualization, we put forward the next step to color the inner hexagons of the 
ruby lattice: a colored link, say red, connect the red inner hexagons. Accordingly, other 
inner hexagons and links can be colored, and eventually we are left with a colored lattice. 
The emergence of the topological color code is beautifully pictured in FigJHl Geometrically, 
it corresponds to shrinking the blue triangles of the original lattice into points, which will 
be referred as sites of a new emerging lattice, see FigJH] (left). Thus, we realize that the 
inner hexagons and vertices of the model are colorable, see FigJH] (middle): if we regard blue 
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(a) (b) (c) 

FIG. 8: (color online) The three stages showing the emergence of the topological color code: (left) 
the original lattice for the 2-body Hamiltonian ([6]). The colors in the links denote the type of spin- 
spin interactions; (middle) a different coloring of the lattice is introduced based on the property 
that the hexagons are 3-colorable, as well as the vertices; (right) the hexagonal lattice obtained 
by shrinking to a point the blue triangles of the original lattice, which become sites in the final 
hexagonal lattice. This corresponds to the strong coupling limit in (jlSp . 

triangles as the sites of a new lattice, we get a honeycomb lattice, see FigJH] (right). In fact, 
the model could be defined for any other 2-colex, not necessarily a hexagonal lattice. 

Connection to the 2-colexes can be further explored by seeing how strings on the ruby 
lattice correspond to the colored strings on the effective honeycomb lattice. To this end, 
consider a typical string- net on the ruby lattice as shown in FigJU^a) . This corresponds to a 
non-perturbative picture of the lOMs of the model. The fat parts of the string-net connect 
two inner hexagons with the same color. In this way, the corresponding string-net on the 
effective lattice can be colored as in Figj9]^b). The color of each part of the string- net on 
the effective honeycomb lattice is determined by seeing which colored inner hexagons on the 
ruby lattice it connects. Three colored strings cross each other at a branching point, and its 
expression in terms of Pauli matrices of sites are given by product of Pauli operators written 
adjacent to the sites. How they are determined, will be clear soon. 

It is possible to use colors to label the closed strings on the honeycomb lattice. Before 
that, let us use a notation for Pauli operators acting on effective spins of honeycomb lattice 
r°, where a = x,y,z. We indicate the labels a as c|c := z,c\c := x,c\c := y, where we are 
using a bar operator. To each c-plaquette, we attach three operators each of one color. Let 
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Bj denotes such operators, where low and up indices stand for c-plaquette / and color of 
the closed string attached to the plaquette, respectively. With these notations, the plaquette 
operators read as follows: 

4=n<''^ (10) 

vef 

where the product runs over all vertices of the c-plaquette / in the honeycomb lattice in 
FigJHl^c). Thus, we can write the explicit expression of operators as follows: 





B) = 


vef 


B) = 


B) = 


vef 


Bf = 


-B) 


= n 

vef 



r:\\ (11) 

All these plaquette operators are constant of motions. Again, We can realize the gauge 
symmetry Z2 x Z2 through the relation B^B^j-Bj = 1. On a compact manifold, for example 
on the torus, all plaquettes are not independent. They are subject to the following constraint: 

^5^ = (-l^/^ (12) 

feA 

where the product runs over all plaquettes / in the lattice A, and is the total number of 
plaquettes. 

We can also realize noncontractible strings on the effective lattice which are rooted in 
the topological degeneracy of the model. They are just the lOMs in FigJ7]when reduced on 
the effective honeycomb lattice. Once the inner hexagons of ruby lattice are colored, they 
correspond to colored strings as in Figj9l Let S*^ stands for such string, where indices fi and 
c denote the homology and color of the string, respectively. This string operator is tensor 
product of Pauli operators of qubits lying on the string. Namely, the string operator is 

V 

The contribution of each qubit is determined by the color of the hexagon that the string 
turns on it, see FiglTOl For example in the string Si shown in FigJTOl the color of the 
plaquettes appearing in f|T3|) marked by light circles. With this definition for string operators, 
the contribution of Pauli operators in the string-net on the effective lattice in Fig|9](b) are 
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(a) 



(b) 



FIG. 9: (color online) An illustraion of correspondence between (a) strings on the ruby lattice, 
corresponding to a non-perturbative picture, and (b)colored strings on the effective honeycomb 
lattice. 

reasonable. Non-contractible colored strings are closely related to the topological degeneracy 
of the model, since they commute with color code Hamiltoniand^]) being integrals of motion, 
but not always with each other. In fact, two strings differing in both homology and color 
anticommute, otherwise they commute. For example let us consider two non-contractible 
closed strings 5*1 and 5*2 corresponding to different homologies of the torus. As shown in 
FigHni they share two qubits. First, suppose both strings are of blue type. The contribution 
of Pauli operators of these two qubits in string 5*1 is , while for string 5*2 the contribution 
is Tj^T"! implying [S\, S2] = 0. Then, let S2 be a green string. In this case the contribution 
of qubits will be rf r|, which explicitly shows that {S^, } = 0. The interplay in ([7]) can be 
translated into an interplay between color and homology as follows. 



where s is the number of sites on the string. This interplay makes the ground state subspace 
of the color code model be a string-net condensed phase. 

B. Degenerate Perturbation Theory: Green Function Formalism 

In this subsection we put the above correspondence between original 2-body lattice Hamil- 
tonian and color code model on a quantitative level. In fact, there is a regime of coupling 
constants in which one of the phases of the 2-body Hamiltonian reproduces the TCC many- 




(14) 
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FIG. 10: (color online) A piece of effective lattice. The strings 5*1 and S2 correspond to different 
homology classes of the manifold. Their expression in terms of Pauli operators are given by their 
associated color and the fact that how they turn on plaquettes on the lattice. 

body structure and physics. In particular, we show that this corresponds to the following 
set of couplings in the original 2-body Hamiltonian 

Jxi Jyi Jz ^ 0, Jxi Jy ^ Jzi (-^^) 

that is, a strong coupling limit in J^. The topological color code effectively emerges in 
this coupling regime. This can be seen using degenerate perturbation theory in the Green 
function formalism. Let H = Hq + V he a. Hamiltonian describing a physical system with 
two-body interaction, and we regard the ||^||, the norm of V", be very small in comparison 
with the spectral gap of unperturbed Hq. We also suppose that Hq has a degenerate ground- 
subspace which is separated from the excited states by a gap A. The effect of V will be 
to break the ground state degeneracy in some order of perturbation. Now the interesting 
question is whether it is possible to construct an effective Hamiltonian, Hes, which describes 
the low energy properties of the perturbed Hamiltonian H. The effective Hamiltonian 
arises at orders of perturbation that break the ground state degeneracy. From the quantum 
information perspective the Hamiltonian H acts on the physical qubits while the effective 
Hamiltonian acts on the logical qubits projected down from the physical qubits. 

We will clarify that how many-body Hamiltonian in ([5]) will present an effective descrip- 
tion of low lying states of the 2-body Hamiltonian ([6]). We use the perturbation about the 
Hamiltonian in ([6]) considering the coupling regim (ITS!) . Here, the qubits on the triangles are 
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physical qubits, and logical qubits are those living at the vertices of the 2-colex. We refer to 
triangles as sites, since they correspond to the vertices of the 2-colex. Thus a triangle will 
be shown by index v and its vertices by Latin indices In fact the low lying spectrum of 
2-body Hamiltonian encodes the following projection from the physical qubits to the logical 
ones at each site: 

n = i^)(ttti + i^)aui, (16) 

where | ft) and | \\) stands for the two states of the logical qubit at site and 1 1) (I i)) is 
usual up (down) states of a single spin in computational bases. 

To this end, we split the 2-body Hamiltonian into two parts. The unperturbed part is 
Hq = — Jz Ylih-iink ^i'^j- limit of strong Ising interaction the system is polarized. The 

interactions between neighboring qubits on different triangles are included in V. They are 
erf (tJ and erf crj corresponding to red and green hnks in Figj3l respectively. So, the transverse 
part of the Hamiltonian is 

r—link g—link 

In the case when ^ Jx, Jy the low lying excitations above the fully polarized state can 
be treated perturbatively. 

The unperturbed part of the Hamiltonian, Hq, has a highly degenerate ground space 
because, for each triangle, two fully polarized states | ttt) and | |4,i) have same energy. The 
ground state subspace is spanned by all configurations of such polarized states. Let N be 
the number of triangles of the lattice. The ground state energy is = -mJz and the 
dimension of the ground space of the Hq or ground state degeneracy reads Qq = 2^ . The first 
excited state is produced by exciting one of triangles and has energy 
with degeneracy gi = QN2^~^. The second excited state has energy Ef^ = (-3A^ + 8) 
with degeneracy g2 = 18N{N — 1)2^"^, and so on and so forth. 

We analyze the effect of V on the ground state manifold by using the degenerate pertur- 
bation theory jsil] in couplings Jx and Jy. We are interested in how ground state degeneracy 
is lifted by including the interaction between triangles perturbatively. Let C stand for the 
ground state manifold with energy E^^ and let P be the projection onto the ground state 
manifold C. The projection is obtained from the degenerate ground states as follows: 

p = llp. , n = i^)(ttti + iu)aui- (18) 
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Using the projection and Green's function we can calculate the effective Hamiltonian at any 
order of Perturbation theory. The eigenvalues of the effective Hamiltonian appear as 
the poles of the Green function G{E) — P[1/{E — H)]P. The effect of perturbation can be 
recast into the self-energy by expressing the Green's function as 1/{E — E^^ — T,{E)). 
So, the effective Hamiltonian will be 

oo 

^eff = E^e2=^i°^+^(^)- (19) 
1=0 

The self-energy can be represented in terms of Feynman diagrams and can be computed for 
any order of perturbation: 

oo 

j:{e) = pyJ^^^'p^ (20) 

n=0 

where U = [1/{E — -f^o)](l ~ P)V. The energy E can also be expanded at different orders 
of perturbation, E = E^^ + X^z^i -^o Now, we are at the position to determine different 
orders of perturbation. Each term of V acts on two neighboring physical qubits of different 
triangles. At a given order of perturbation theory, there are terms which are product of 
(7^ and (T^ acting on the ground state subspace. Each term when acts on the ground state 
manifold brings the ground state into an excited state. However, there may be a specific 
product of the and which takes the ground state into itself, i.e. preserve the polarized 
configurations of triangles. 

At zeroth-order the effective Hamiltonian will be trivial H^^ — E^K The first-order 
correction is given by the operator 

H^^ = PVP (21) 

The effect of V is to move the states out of the ground state manifold because each term either 
(jx^x (j2/(j2/ flip two qubits giving rise to two triangles being excited, i.e VP — P2VP, where 
the operator P2 is the projection to second excited state manifold. Therefore, PVP — 0, 
and there is no first-order correction to the ground state energy. 

The second-order correction to the ground state will be the eigenvalues of the following 
operator. 

H^l^ = PVG'{Ef^)VP + PVP, (22) 
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where the operator G'{E) = 1/{E — Hq) is the unperturbed Green's function and the super- 
script prime stands for the fact that its value be zero when acts on the ground state. The 
second-order correction only shifts the ground state energy, and therefore, the second-order 
effective Hamiltonian acts trivially on ground state manifold, 

^i? = 3^4^^- (23) 

In fact the first V flips the qubits and the second V flips them back. As we go to higher 
order of perturbation theory the terms become more and more complicated. However, if 
the first-order is zero as in our case, the terms becomes simpler. Thus, the third-order of 
perturbation will be zero and will leave corrections to energy and ground state intact: 

H^^ = PV (G'{Ef^)V) ' P = 0. (24) 

The forth-order of perturbation theory contributes the following expression to the correction 
of ground state manifold: 

h'^ = PV (G'{ES'^)Vy P - Ef^PV (G\Ef^)yVP, 

(25) 

where Eq is the second order correction to the ground state energy obtained in (l23|l . The 
first term includes four V and must act in the ground state in which the last V returns 
the state to the ground state manifold. The second term is like the second-order. There 
are many terms which must be calculated. However, since the forth-order only gives a 
shift to the ground state energy, we don't need them explicitly. So, we can skip the forth- 
order. Fifth-order correction yields terms each containing odd number of V , so it gives zero 
contribution to the effective Hamiltonian. 

The sixth-order of perturbation leads to the following long expression. 

H^^ = PV (G"(Ei°))r)' P - Ef'PV (G'iES'^)yVP 
+ (Ef)'py (G'(Er))Vp 
-Ef PV (G'(Ei°)))V (G"(Ef )v)'p 
-EPpVG'{eS''^)V (G\Ef'^)yVG'{Ef^)VP 
-Ef^PV (g"(eJ°V)' (G"(^r))' (26) 
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Apart from the first term, other terms contain two or four V and as we discussed in the 
preceding paragraphs they only contribute a shift in the ground state energy. However, the 
first term gives the first non-trivial term breaking in part the ground state degeneracy. In 
the sixth order correction, there are some terms which are the product of a^a^ and cr^cr^ 
associated to the red and green links of the ruby lattice. Some particular terms, as seen 
below, may map ground state subspace into itself. For instance, consider the following 
product of links around an inner hexagon 

n « = ± n (27) 

<i,j> ieVft 

where the first product runs over three red and three green links making an inner hexagon, 
Vh stands for the set of its vertices and the prefactor ± depends on the ordering of links in 
the product. The action of a on one vertex (qubit) of a triangle encodes an logical 
operator acting on the associated vertex of lattice A. This can explicitly be seen from the 
following relation: 

< = P,a^P, = |^)(^|-|^)(U|, (28) 

where acts on one of the vertices of a triangle and P„ is the projection defined in ( flGjl . 
Thus, the expression of (1271) can be related to the plaquette operator = — H'^i?' where 
the index / denotes a plaquette of effective lattice A as in FigJT] and product runs over six 
sites around it. Now we go on to pick up the sixth-order correction to the ground state 
manifold. There are many terms which must be summed. Sixth-order correction up to a 
numerical constant contributes the following expression to the effective Hamiltonian: 

= constant - 6^ (29) 

f 

where 5 is a positive numerical constant arising from summing up 720 terms related to 
the order of product of six links around an inner hexagon of ruby lattice. Although, its 
exact numerical value is not important, but knowing its sign is essential for our subsequent 
discussions. As it is clear from the first term in fj26|) . five Green's functions and six V in 
the perturbation contribute a minus sign to the expression. This minus sign together with 
the sign appearing in (^7^ enforce the coefficient 5 be a positive constant. Now it is simple 
to realize how the vectors in the ground state manifold rearranged. Trivially, all plaquette 
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operators commute with each other and their eigenvalues are ±1. All polarized vectors 
in C are the eigenvector of the effective Hamiltonian emerging at sixth-order. But, those 
are ground states of the effective Hamiltonian in ( l29l) which are eigenvectors of all plaquette 
operators with eigenvalue +1. Thus, highly degenerate ground state of the unperturbed 
Hamiltonian is broken in part. The same plaquette operators Bj also appear at higher 
order of perturbation. For example, at eighth-order. Instead of giving the rather lengthy 
expression of the eighth-order correction, we only keep terms resulting in the plaquette 
operators as follows: 

t5 73 I 73 75 

HjS = constant - /3 ^ ^ ^ Yl (^0) 

where /3 > 0. This term is added to the one in (|29l) to give the effective Hamiltonian up 
to eighth-order, but the ground state structure remains unchanged. Further splitting in the 
ground state manifold is achieved by taking into account the ninth-order of perturbation. 
The expression of ninth-order is very lengthy. However, the first term of the expression 
containing nine V gives some terms being able of mapping the ground state manifold into 
itself in a nontrivial way. These terms map a polarized triangle, say up, to a down one. 
Indeed, when one or two qubits of the polarized triangle gets flipped, its state is excited. 
However, flipping three qubits of the triangle returns back the ground state onto itself. This 
process encodes and logical operators acting on logical qubits arising through the 
projection. Let af, af and af act on three qubits of a triangle. The encoded operator 
will be 

ry = P^ala^alP, = ^)(^ | + ^| ^) |. (31) 

If a^, cr| and (j| act on three qubits of a triangle, the encoded logical operator will be 

- r: = P.a^faf = _ | ^) | _ | ^) | . (32) 

As we already pointed out a plaquette of the ruby lattice is made up of an inner hexagon, 
an outer hexagon and the six blue triangles. It is possible to act on the polarized space of 
the blue triangles by making two different combinations of 9 link interactions: i/ Applying 
6 link interactions on the outer hexagon (three of them of XX type and another three of YY 
type), times 3 link interactions of red type on the inner hexagon. Notice that every vertex of 
the blue triangles in the plaquette gets acted upon these 9 link interactions. The resulting 
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effective operator is of type due to (|3T|) . ii/ Applying 6 link interactions on the outer 
hexagon, times 3 link interactions of green type. Then, the resulting effective operator is of 
type due to ([321). 

The effective Hamiltonian at this order then reads 

= PV (G"(4°y)' P + ... (33) 

= constant - 7^ E ^1 " E ^f- 

Again, the sign of coefficient 7 is important. Nine V^s, six or r^, and eight Green's 
function imply that the 7 must have positive sign. 

Putting together all above obtained corrections lead to an effective Hamiltonian encod- 
ing color code as its ground state[l|, Q. Therefore, up to constant terms, the effective 
Hamiltonian reads as follows 

/ / / 

where the kz, kx and ky are positive coefficients arising at different orders. Since BJBj = Bj, 
the above effective Hamiltonian is just the many-body Hamiltonian of the color code as in 
([1]). The terms appearing in the Hamiltonian mutually commute, so the ground state will 
be the common eigenvector of plaquette operators. Since each plaquette operator squares 
identity, the ground state subspace, C, spanned by vectors which are common eigenvectors 
of all plaquette operators with eigenvalue +1, i.e 

C = m: B'j\^|J) = lij), Bjl^) = \^) ; V/}. (35) 

The group of commuting boundary closed string operators can be used as an alternative 



way to find the terms appearing in the effective Hamiltonian [82]. As we pointed out in 
the preceding section, the non-zero contribution from various orders of perturbation theory 
results from the product of red and green links which preserve the configurations of the 
polarized triangles, i.e maps ground state manifold onto itself. For instance consider the 
elementary plaquette operator A corresponding to a closed string in FiglSl Each triangle 
contributes a^a^a^ to the expression of operator which is projected to as in (!32|) . Thus, 
the effective representation of plaquette operator reads as follows: 

PSfP B]. (36) 
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Plaquette operators Sj and 5*^ can also be recast into the effective forms as follows: 

PSfP Bj , PSfP ->-B}. (37) 

These are lowest order contributions to the effective Hamiltonian as we obtained in (IMll . 
Higher order of perturbation will be just the product of effective plaquette operators. The 
nontrivial strings winding around the torus will have also effective representations and appear 
at higher orders of perturbation. In general, every string-net lOM on the ruby lattice are 
projected on an effective one as in Figj9l 



V. BOSONIC MAPPING 

As we stated in Sect JIIIl one of the defining properties of our model is the existence of 
non-trivial integrals of motion lOM, called string-nets. As a particular example, the Kitaev's 
model on the honeycomb lattice has strings lOMs, but not string-nets. We are interested in 
models in which the number of lOMs is proportional to the number of qubits in the lattice, 
i.e., 

/ = vNg, (38) 
where / is the number of lOMs, Ng is the number of qubits (spins) in the lattice, and rj 



is a fraction: = ^ for the Kitaev model 72|, ''7 = | for our model in ([6]) [lO(]. The fact 



that ?7 is a fraction t] < 1 implies that these models based on string-net lOMs will not 
necessarily be exactly solvable. In the Kitaev's model, it turns out to be exactly solvable 
using an additional mapping with Majorana fermions, but this need not be a generic case. 
Therefore, we need to resort to other techniques in order to study the physical properties of 
these models. We consider here and next section, approximate methods based on bosonic 
and fermionic mappings. The application of the bosonic method to our model is based on a 
mapping from the original spins on the blue triangles to hardcore bosons with spin [l^. With 
this mapping, it is possible to use the PCUTs approach (Perturbative Continuous Unitary 



Transformations) 83(|. This is inspired by the RG method based on unitary transformation 



introduced by Wegner (the Similarity RG method) 8J]. Originally, the PCUTs method was 



applied to the Kitaev model [85[. As we will see below, the method paves the way to go 
beyond the perturbation approach presented in the preceding section which fits into a sector 
without any hard-core boson. The physics at other sectors is very promising and we study 
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it here. Emergence of three famihes of strongly interacting anyonic fermions and invisible 
charges are among them, which are not present in the Kitaev's model or its variants. To start 
with, let us set for simplicity = 1/4, and consider the extreme case of Jx = Jy = 0. In 
this case the system consists of isolated triangles. The ground states of an isolated triangle 
are polarized states | ttt) and | ^ii) with energy —3/4. The excited states that appear by 
flipping spins are degenerate with energy 1/4. In this limit, the spectrum of whole system 
is made of equidistance levels being well-suited for perturbative analysis of the spectrum: 
Green's function formalism as discussed in the preceding section, which may capture only 
the lowest orders of perturbation or another alternative approach based on the PCUT. The 
change from the ground state to an exited state can be interpreted as a creation of particles 
with energy +1. This suggest an exact mapping from the original spin degrees of freedom 
to quasiparticles attached to effective spins. The mapping is exact, i.e. we don't miss any 
degrees of freedom. Such a particle is a hard-core boson. At each site, we attach such a 
boson and also an effective spin-|. Let choose the following bases for the new degrees of 
freedom 

\a,d) = \a) ® \d), a =ff, J|, ci = 0,r, g, b, (39) 

where a and d stand for states of the effective spin and quasiparticle attached to it. The 
Hilbert space Tic representing the hard-core bosons is four dimensional spanned by bases 
{|0), |r), |g), |b)}. Now the following construction relates the original spin degrees of freedom 
and new ones in (!39l) 



^,0) = 


Ittt), 


1^,0) = 


liii) 




im), 




iin) 




liti), 




Itit) 




lUt), 




im) 



(40) 



Within such mapping, the effective spins and hard-core bosons live at the sites of the effec- 
tive hexagonal lattice A in Fig. [H](c). Recall that this lattice is produced by shrinking the 
triangles. At each site we can introduce the color annihilation operator as be '■= |0)(c|. The 
number operator n and color number operator Uc are 

n:=^nc, := blbc- (41) 
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Annihilation and creation operators anticommute on a single site, and commute at different 
sites, that is why they are hard-core bosons. We can also label the Pauli operators of original 
spins regarding to the their color in Fig. Mjo) as with c = r, g, b. The mapping in ( HOll 
can be expressed in operator form as follows 

= ® p„ < = r'^ ® (6j + bc + s^r,) , (42) 

where i' = x,y, := —Sy := 1, the symbols r denote the Pauli operators on the effective 
spin and we are using the color parity operators Pc and the color switching operators Vc 
defined as 

: = 1 - 2 (ng + ng) , : = blbg + bibs- (43) 

Now we can forget the original ruby lattice and work on the effective lattice in which the 
bosons are living at its sites. With the above identification for Pauli operators, the 2-body 
Hamiltonian in (jH]) can be written in this language. Before that, let fix a simplified notation. 
All spin and bosonic operators act on the sites of the effective lattice. We refer to a site 
by considering its position relative to a reference site: the notation O c means O applied at 
the site that is connected to a site of reference by a c-link. The 2-body Hamiltonian then 
becomes 

H = -3N/4 + Q - ^ ^ J,,|, T/, (44) 

A c^c' 

with the number of sites, Q := "^^n the total number of hardcore bosons, the first 
sum running over the A^ sites of the reduced lattice, the second sum running over the 6 
combinations of different colors c, c' and 

=ui + ^^ + f + h.o., (45) 

a sum of several terms for an implicit reference site, according to the notation convention 
we are using. The meaning of the different terms in fH5|) is the following. The operator 
is a c-boson hopping, switches the color of two c- or c-bosons, -u^' fuses a c-boson with a 
c-boson (or a c-boson) to give a c-boson (c-boson) and destroys a pair of c-bosons. The 
explicit expressions are 

ti ■ = ' bcblc' ' ■= ' rcTcc' , 

:= Sc'\cTcbcrc,c', vi := r^'bcb^c', (46) 
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where we are using the notation 

ri := r^'lv;;'^ (47) 

We can also describe the plaquette lOM operators in Fig. [5] in terms of spin-boson degrees 
of freedom by means of the mapping in fl42l) . For each plaquette / and color c, the plaquette 
operator is expressed as 

S'f--=\{rfPc'.c. (48) 

where c' is the color of the plaquette /, the product runs through its sites and * is just a 
convenient symmetric color operator defined by 

CT*rc:=c, c-k c := c-k c := c. (49) 

The relation in fllS]) is just a generalization of plaquette operators in flTT]) to other sectors 
of the system. In fact taking the zero particle sector, the expressions in (fTTl) are recovered. 
In the same way the nontrivial string operators in FigJT] can be described with the above 
mapping as 

S'',--=\{rfpc.c'. (50) 

where denotes the homology class of the string. On closed surfaces, not all plaquette 
operators are independent. They are subject to the following constraints 

^^^ = (-l)^/^ n ^/ = (-lr/^ (51) 

/eA c=r,g,b 

where s is the number of sites of a given plaquette /. The first equation can be further 
divided into products over subsets of plaquettes giving rise to other constants of motion, the 
so called color charg 

n^/n^j= n^/n^/=np^ (52) 

f&f^a /eAg /eAg feKa A 
In these products the spin degrees of freedom are washed out, since they appear twice and 
consequently square identity. By use of equation ( 143|) . the product over parity operators can 
be written as 

XA(c) = ^Pc = (-l)^^"+'^^ (53) 

A 

where Qc = ^a'^c is the total number of c-bosons. It is simple to check that the above 
equation commutes with Hamiltonian in (jH]). For each family of bosons we can attach a 
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charge. We suppose that each c-boson carries a charge as Xc that is, an irrep of the gauge 
group. In particular, the Hamiltonian preserves the following total charge 

XA ■■= xf'X^'X^'- (54) 
A. Emerging particles: anyonic fermions 

Equation fl52l) could already suggest that the parity of vortices are correlated to the parity 
of the number of bosons. In particular, creation of a c-boson changes the vorticity content 
of the model. 

The statistics of vortices depend on their color and type as in FigJS But what about the 



statistics of c-bosons? As studied in [SG], the statistics of quasiparticles can be examined 
using the hopping terms. These hopping terms are combined so that two quasiparticles 
are exchanged. In addition to usual hopping terms, we also need for composite hopping , 
that is, a c-boson hops from a c-plaquette to another, which is carried out by terms like 
f: = u^up,. = u'^T.ui\. Let us consider a state with two c-boson excitations located at two 
different sites separated from a reference site by c and c links. An illustrative example for the 
case of, say, blue bosons is depicted in FiglTTTa). Consider a process with the net effect of 
resulting into the exchange of two bosons. Each step of process can be described by hopping 
terms. Upon the combination of hopping terms, we are left with the following phase 

iKA^^^^^ = ir'ryrXr:^)' = -1, (55) 
which explicitly show that the quasiparticles made of hardcore-bosons and effective spins 



have fermion statistics [lOl]. Thus we have three families of fermions each of one color. 
These are high energy fermions interact strongly with each other due to the fusion term in 
the Hamiltonian. Fermions from different families have mutual semionic statistics, that is, 
encircling one c-fermion around a c-fermion picks up a minus sign. This can also be checked 
by examining the hopping terms as in FigJTlTa). Thus we are not only dealing with fermions 
but also with anyons. 

The elementary operators in (H6l) have a remarkable property, and that, they all commute 
will plaquettes and strings lOM. This naturally implies any fermionic process leaves the 
vorticity content of the model unaffected. A fermionic process may correspond to hopping, 
splitting, fusion and annihilation driven by the terms in the Hamiltonian (145!) and P6|) . A 
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FIG. 11: (color online) (a) A step by step illustration of hopping of two blue fermions with the net 
result of exchange of two fermions (b)A representation of a fermionic process in which a c-fermion 
(here red) at the origin site is annihilated and undergoes hopping, splitting and fusion processes, 
and then is created at another site. 



typical fermionic process is shown in Fig JTTT b). with the net result of displacement of a 
r-fermion from a site at the origin to other one of the effective lattice. A very feature of 
this process is the existence of vertexes, which is essential to high-energy fermions, that is, 
three fermions with different color charge can fuse into the vacuum sector. At the vertex, 
three different colored strings meet. Notice that the colored strings shown here have nothing 
to do with the ones we introduced in Sect JIIII Indeed, these are just product of some red 
and green links of the ruby lattice. When translated into the spin-boson language, they are 
responsible for transportation of c-fermions through the lattice. 

Now we can think of constraints in (l52l) and ( l53l) as a correlation between the low energy 
and high energy sectors of our model. They explicitly imply that the creation of a c-fermion 
creates the vortices with net topological charge of (c, Xs)- Alternatively, as suggested by 
mapping in ( HOj) . flipping of spin on a triangular can create or destroy the excitation, that 
is, a high energy fermion can be locally transformed into low energy ones. This amounts 
to attach a topological charge from low energy sector to high energy excitations. Thus a 
c-fermion carries a topological charge. On the other hand, an open c-string commutes with 
all plaquettes except some of them, so they create or destroy a particular charge among the 
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charges in Figj2l It is simple to check that which charge they carry at their open ends. In fact 
they carry (c, Xc) charges. These latter charges have trivial mutual statistics relative to the 
charges carried by high energy fermions, since they belong to different families of fermions. 
Thus the charge carried by a c-string must be invisible to the high energy fermions. 

B. Perturbative continuous unitary transformation 

The physics behind Hamiltonian in flU]) can be further explored by resorting to approxi- 
mate methods. The specific form of the energy levels of our model in the isolated limit, the 
existence of equidistant levels, makes it suitable for perturbative continuous unitary trans- 
formations. In this method the Hamiltonian is replaced by an effective one within unitary 
transformations in which the resulting effective Hamiltonian preserves the total charges, i.e. 
[H'^^,Q] = 0. Thus the analysis of the model relies on finding the effective Hamiltonian 
in a sector characterized by the number of charges at every order of perturbation. For our 
model each sector are determined by the number of c-fermions. Each term of the effective 
Hamiltonian in any sector is just a suitable combination of expressions in fH6|) in such a way 
that respects the total color charge in that sector. For now we briefly analyze the lowest 
charge sectors. 

In the zero-charge sector, only the effective spin degrees of freedom do matter. The 
effective Hamiltonian is just a many-body Hamiltonian with terms that are product of 
plaquette operators as follows 

{c} {/} 

where the first and second sum run over an arbitrary collection of colors and plaquettes 
of effective honeycomb lattice. The coefficients O's are determined at a given order of 
perturbation. The product of plaquette operators is nothing but the string-net operators. 
Let us focus at lowest order of perturbation, where the model represents non-interacting 
vortices. First, let us redefine plaquette operators as 

Bf=jfll<^ B'f=jfll^l (57) 

where := Jw/\Jw\- At ninth order of perturbation the effective Hamiltonian is 

^0*" = - E i^-^f + ^y^} + kzB]Bj) , (58) 
/eA 
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with [10[ 

...2|,,,/ + 0(n ^ . . 5|||,.,,|3 (59) 

This is exactly the many-body Hamiltonian of topological color code obtained in flM|) using 
degenerate perturbation theory with the additional advantage of knowing the coefficients 
exactly. Its ground state is vortex free and can be written explicitly by choosing a reference 
state as 

2^/^-^n(^)l^)®|0V (60) 

Other degenerate ground states can be constructed by considering the nontrivial string 
operators winding around the torus. Excitations above ground state don't interact. Going 
to higher order of perturbation, as equation in fl56|) suggests, the ground states remain 
unchanged, however the excitation spectrum changes and vortices interact with each other. 

The one-quasiparticle sector can also be treated by examining the expressions in ( l46l) . 
The effective Hamiltonian can be written as 

Hf = Hf-J20RRblRbc. (61) 

{R} 

What the second term describes is nothing but the annihilation of a c-fermion at a reference 
site and then its creation at a site connected to the reference by a string-net R, as shown 
in Fig JTTT b). Again notice that this string-net is just the product of green and red links 
of original 2-body Hamiltonian, in its effective form is given in terms of spin-boson degrees 
of freedom. The coefficients O's are determined at any order of perturbation. Notice that 
these coefficients are different from those in (|56l) . In the first order, only the hopping term 
does matter. Let us consider the sector containing a c-fermion. Up to this order, the fermion 
can only hop around a c-plaquette. This implies that at first order the fermion perform an 
orbital motion around a plaquette of its color. Notice that the fermion can not hop from 
a c-plaquette to other c-plaquette at the first order, since it needs for a composite process 
which appears at second order. This composite process is a combination of splitting and 
fusion processes. This is a virtual process in the sense that the splitting of a c-fermion into 
two c- and c-fermion takes the model from 1-quasiparticle sector into the 2-quasiparticle 
sector, but the subsequent process fuses two particles into a single one turning back to 1- 
quasiparticle sector. Thus, at second order the c-fermion can jump from one orbit to other 
one. 
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At first order, for J = = Jy, we get a — 2J contribution to tfie energy gap coming 
from orbital motion. Going to second order we get a non-fiat dispersion relation. The gap, 
at this order, is given by 1 — 2J — and thus it closes at J ~ 0.45. This is just an 
approximate estimation, since we are omitting all fermion interactions and, perhaps more 
importantly, we are taking J ~ J^. However, it is to be expected that as the couplings 
Jx ~ Jy grow in magnitude the gap for high-energy fermions will reduce, producing a phase 
transition when the gap closes. Such a phase transition resembles the anyon condensations 



discussed in 



87H89|]. There are three topological charges invisible to the condensed anyons. 



This means that in the new phase there exists a residual topological order related to these 
charges. They have semionic mutual statistics underlying the topological degeneracy in the 
new phase. 



C. Fermions and gauge fields 



The emerging high-energy c-fermions always appear with some nontrivial gauge fields 10|, 



861], and carry different representation of the gauge symmetry Z2 x Z2 of the model. Before 
clarifying this, we can see that the plaquette degrees of freedom correspond to Z2 x Z2 gauge 
fields. This correspondence is established via introducing the following plaquette operators 

d5 As/2 qc 

3x ^f, 

B} := H.jyy/' S}. (62) 

The gauge element g/ G Z2 x Z2 that can be attached to the plaquette / is determined by 
following eigenvalue conditions 

Xc(g/) = B'}, (63) 

which always has a solution due to 

{B'}Y = B)Bp) = 1. (64) 

The ground state of color code Hamiltonian ( [58|) is vortex free and corresponds to Xc = 1- 
The fact that for a 2-colex with hexagonal plaquettes, the gauge fields can be related to 
representation of the group is immediate. One way to see this is to check the phase picked 
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FIG. 12: (color online) A piece of 2-body lattice corresponding to a 2-colex composed of octagon 
and square plaquettes. The lattice is 3-colorable as in FiglH) 

up by a c-fermion when it moves around a plaquette. Turning on a plaquette is done by 
combination of hopping operators yielding the phases as Bj , Bj and Bj that are consistent 
with ( l64l) . However, this is not generic for all 2-colexes. For 2-colexe plaquettes that the 
number of their edges is a multiple of four, we see that the ground state carries fluxes. 
Perhaps the most important of such lattices is the so called 4-8-8 lattice shown in FigJT2l It 
contains inner octagons and squares. Once degenerate perturbation theory is applied about 
the strong limit of the system, the effective color code Hamiltonian in terms of plaquette 
operators in fl62|) at 12th order of perturbation is produced, as follows 

Heff = - '^{k.^Bj + kyBj + k^Bj) + multiplaquette terms, (65) 
/ 

where sum runs over all squares and octagons. Notice that at 12th order of perturbation 
multiplaquette terms that are product of square plaquette operators are also appeared. It is 
simple to check that the coefficients fc's have positive sign. As we can relate the plaquettes 
to the representation of gauge group, the ground state corresponds to vortex free sector. 
In fact, the ground state of all 2-colexes with plaquettes of any shape is vortex free and 
correspond to Xc = 1 of gauge group. What is able to differentiate between ground states 
of 2-colex plaquttes with 4n (n an integer) edges from others is related to the gauge fields 
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FIG. 13: (color online) A c-fermion process. A red fermion goes around a region A. The phase 
that it picks up depends on vortex configurations shown by spiral lines, blue and green fermions 
and number of plaquettes that the number of their edges is of multiple four. 



attached to a fermion. In particular, there is a background of 7r-fluxes in the ground states of 
such lattices, and the emerging c-fermions can detect them. To make sense of the existence 
of such fluxes, let us consider a simple fermionic process as explained above. When a c- 
fermion turns on a plaquette, the combination of hopping terms yield —Bj, —Bj and —Bj, 
which clearly imply that 



-Sj)(-i?f)(-i?^) = -l. 



(66) 



This result exhibit that the ground states of color code models defined on lattices with 4n- 
plaquettes carry fluxes. Thus, in its representation in terms of Z2 x Z2 gauge group, the 
fluxes must be subtracted away. 

Our derivation for ground states with 0- and vr-flux can be compared with Lieb 
theorem tgci]. which states for a square lattice the energy is minimized by putting tt flux 
in each square face of the lattice. The connection to our models makes sense when we con- 
sider how 2-colexes with hexagonal and 4-8-8 plaquettes can be constructed from a square 
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lattice by removing some edges. Then the total vr fluxes in a set of square faces correspond- 
ing to a 2-colex plaquette amounts to the flux that it carries. It is simple to see that each 
hexagonal plaquette composed of two (imaginary) square faces, and flux vr in each square 
face then implies flux in the hexagon. The same strategy holds for fluxes carried by the 
4-8-8 plaquette(and in general for all 4n-plaquettes) . Once again we see that each plaquettes 
of latter 2-colxes is composed of odd number of square faces, thus they carry flux vr in their 
plaquettes. 

Now we can give a general expression for the gauge fields seen by emerging high-energy 
fermions. To do so, let us consider a process in which a c-fermion is carried around a region 
A, as in FigjlSl The hopping process yields a phase 



where q\ = Yl^ qj, denotes the number of c-fermions inside A and the number of 2- 
colex plaquettes inside A with a number of edges that is a multiple of four. Thus, we can see 
that each family of fermions carries a different representation of gauge group given by values 
of inside the region. Moreover, it emphasizes that fermions with different color charges 
have mutual semionic statistics. Clearly for hexagonal lattices = 0, and the ground state 
carries no fluxes. 

VI. FERMIONIC MAPPING 

In this section we will come back to the original Hamiltonian of ([6]) on the lattice in order 
to use another approximate method based on fermionic mappings. This Hamiltonian can 
be fermionized by Jordan- Wigner transformation [73]. To do so, firstly it is convenient to 
present a lattice which is topologically equivalent to the lattice of Fig|3l This is a new type 
of "brick- wall" lattice as shown in FigJTH The black and white sites are chosen such that, 
at the effective level, the lattice be a bipartite lattice, since the effective spins are located at 
the vertices of hexagonal lattice representing a bipartite lattice. Note that neither original 
lattice in Figj3] nor the brick-wall one in FigJH] are bipartite in their own. Also as we will 
see, the fermionization of the model needs some ordering of sites in the "brick- wall" lattice. 
The unit cell of the brick-wall lattice is comprised of two triangles as shown in FigJH] in a 
yellow ellipse. The translation vectors Ui and n2 connect different unit cells of the lattice. 




(67) 
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The deformation of the original lattice into a "brick- wall" lattice allows one to perform 
the one dimensional Jordan- Wigner transformation. The one dimensional structure of the 
lattice is considered as an array of sites on a contour as shown in FigJTSl The sites on the 
contour can be labeled by a single index and the ordering of the sites is identified by the 
direction of the arrows in FigJTSl The expression of the Pauli operators in terms of spinless 
fermions will be: 



'^t ~ 6xp j iir'^^ajai j , crj = {^Ojaj — 1), (68) 

Kj 

where spinless fermions satisfy the usual anticommutation relations as follows 

{tti, a]} = Sij, {al, a]} = {a^, aj} = 0. (69) 
Next, we introduce Majorana fermions as follows: 

Cj = —i{aj — CLj), dj = a] + aj (70) 

for black sites and 

Cj = Oj + ttj, dj = —iia] — dj) (71) 

for white sites. Majorana operators are Hermitian and satisfy the following relations: 

kj = 1, kjki = kikj i ^ j, ('^2) 

where k = c,d. For each unit cell we can realize two vertical red and green links correspond- 
ing to a^a^ and a^a'^ interactions. Observe that each vertex of the lattice can be specified 
by three indices as follows: (i) index I is introduced to specify which unite cell the vertex 
belongs to, (ii) b (w) denotes the black (white) vertex, (iii) /, j or k label the position of the 
vertex in each triangle. With this labeling of vertices, the expression of the different terms 
appearing in the two-body Hamiltonian in terms of Majorana fermions is listed below: 

blue links : crja^ = iUjkdjdk, 
horizontal r — links : crjcr^ = idjdk, 
horizontal g — links : crjcr^ = —idjd^, 
vertical r - links : af^^Vf^^^ = tdf^^df^^F,, 

vertical g - links : aj'^^Vj'^^^ = -id^^^^f^ Fi, (73) 
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FIG. 14: (color online) Brick- wall lattice that is topologocally equivalent to the original lattice 
in FiglSj Triangles with black and white vertices are chosen such that at the logical level they 
correspond to a bipartite lattice. A unit cell of the lattice has been shown by yellow region and 
the whole lattice can be constructed by translation vectors ni and n2. The chains supporting the 
field operator F have been shown by brown ribbons. 



.« -o «- 
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FIG. 15: (color online) Ordering of sites in the brick-wall lattice of Figll4[ Such ordering makes it 
possible to perform a one-dimensional Jordan- Wigner transformation. 



where Ujk = —icjCk and the operator Fi is a non-local operator. Interestingly enough, this 
non-local term has the following expression: 

k-l 

Fi = exp(i7r ^ a]ai) = J]^ icidi, (74) 

l=j l€R 

where we have used the ordering of the brick-wall lattice and R stands for a set of spins 
crossed by two brown ribbons shown in FigHH Observe that for each unit cell we can 
realize such ribbons. The non-local operator can be written as a product of some plaquette 
operators each corresponding to a rectangle (inner hexagon) wrapped by the ribbons. Notice 
that an inner hexagon of the ruby lattice in Figj3] looks like a rectangle in the brick-wall 
lattice. Thus, two ribbons are nothing but a combination of rectangles as in FigJT4l This 
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plaquette operator precisely corresponds to operator Sj^ in FigJS] Let /i,/2,... stand for 
those rectangles (inner hexagons). Thus, we have: 



Fi = l[cidi = Sf^^Sf,^®S[ 



c 



(75) 



l&R 



After all these transformations, we arrive at an exact fermionized Hamiltonian given by 



b—link 



r—link 



g—link 



(76) 



I I 
It is simple to check the following commutation relations: 



H,F^ 



H,U. 



, ^jk 



Fi, Fy 



0. 



(77) 



Since Fi = 1, the nonlocal operators can be replaced by its eigenvalues Fi = ±1 . Thus 
the nonlocal terms appearing in the Hamiltonian which are related to the vertical links of 
the brick-wall lattice can be reduced to quadratic terms. However, the first sum in ( !76|) 
which is related to the triangles cannot be reduced to the quadratic term. This is because 
local fields Ujk corresponding to three links of a triangle anticommute with each other as well 
as with some Fis. Due to these anticommutation, all fields cannot be fixed independently. 
This fact is in sharp contrast with Kitaev's model and its variants. These latter models are 
defined on trivalent lattices. The different fields live at spatially separated links allowing 



a free fermion exact solution 
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73|. The obtained quadratic Hamiltonian describes free 



majorana fermions in a background of Z2 charges. Instead, the lattice of our model is 4- 
valent, a sharp difference that prevents complete solvability and gives rise to very interesting 
features not present in the mentioned models. Note that the fields Ujk are highly interacting 
since on a triangle three fields go to vacuum in the sense that (jjkUkiUij = —1- This is 
resemblance of vertex interaction in high energy fermions that we have seen in Sect JV] with 
the bosonic mapping. However, this latter relation doesn't coincide with the symmetry of 
the model as they do not commute with each other. This will be considered next. 

Thus far, we have considered fields that are present in the Hamiltonian. In what follows, 
we introduce another set of fields which have the Z2 x Z2 symmetry commuting with each 
other and with the Hamiltonian. To this end, consider a plaquette /. As before, by a 
plaquette we mean an outer and an inner hexagon with six triangles between them. Let V/ 
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and Vh stand for sets of vertices of plaquette and inner hexagon, respectively. It is natural 
that Vh C Vf. To each plaquette we attach the three following fields: 

0} = n u n = n ^^^i' (^s) 

j&Vf jef\h v£Vh jeVh 

where by f\h we simply mean Vf — Vh- Each (pf squares identity. They commute with each 
other and with Hamiltonian and Fi: 

= [H, 0?] = (79) 

Also, the fields 0/ are responsible of the Z2 x Z2 gauge symmetry since 0]- ® 0j (g) 0^ = — 1. 
The above fields are related to the plaquette operators. Using the transformations we have 
introduced in ( !73|) . we can fermionize the conserved plaquette operators obtained in FigJS] 
They are associated to the above constants of motion as follows 

Sf = (P}F,, Sf = ct>)F,, Sf = <P}. (80) 

Although the above gauge fields make it possible to divide the Hilbert space into sectors 
in which be eigenspaces of gauge fields (or eigenspaces of plaquette operators), they do not 
allow us to reduce the Hamiltonian in fl76|) into a quadratic form. The Fis can be fixed 
as they commute with the Hamiltonian. But, we are not able to reduce the four-body 
interaction terms in the Hamiltonian into quadratic form. In fact, the anticommutation of 
UjkS on a blue triangle prevents them to be fixed consistently with gauge fixing. 

VII. CONCLUSIONS 

We have introduced a two-body spin-1/2 model in a ruby lattice, see FigJHl The model 
exhibits an exact topological degeneracy in all coupling regimes. The connection to the 
topological color codes can be discussed on the non-perturbative level as well as confirmed 
by perturbative methods. In the former case, on the ruby lattice we realized plaquette 
operators with local Z2 x Z2 symmetry of the color codes. All plaquette operators commute 
with the Hamiltonian and they correspond to integrals of motion. The plaquettes can be 
extended to more complex objects that can be considered as string-nets: non-trivial strings 
with branching points. The nontrivial strings corresponding to the various homology classes 
of the manifold determine the exact degeneracy of the model. For the case of periodic 
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boundary conditions, i.e on a torus, and for each non-contractible cycle of the torus, we 
can identify three nontrivial closed strings. Once each of them is colored, the plaquettes of 
the lattice can be correspondingly colored as in FigJHl For each homology class, they are 
related to each other by the gauge symmetry of the model. The crucial property of these 
strings is that they commute with Hamiltonian but not always with each other. This is 
independent of the regimes of coupling constants of the model. Being anticommuting closed 
nontrivial strings, the model has exact topological degeneracy. To clarify this observation, 
we use perturbation theory to investigate a regime of coupling corresponding to a strong 
coupling limit (triangular limit). In this limit the topological color code will be the effective 
description of the model. The effective representation of the closed loop operators determine 
the terms appearing in the effective Hamiltonian at different orders. 

Unlike the Kitaev's model or any its variants, our model is not integrable in terms of 
mapping to Majorana fermions, to the best of our knowledge. This model is a four-valent 
lattice and gauge fields not always commute with each other. However, we have emphasized 
in Sect illll that the existence of exact integrals of motion (lOMs) at a non-perturbative level 
is far more enriching than demanding exact-solvability of a model. In fact, if the number 
of lOMs is large enough, the model can turn out to be solvable. Thus, fixing plaquette 
operators can not give rise to fix all gauge fields. 

The description of our model in terms of hard core bosons yields very fruitful and interest- 
ing physics of the model. Using a bosonic mapping, it is possible to discuss the emergence of 
strongly interacting anyonic fermions. They form three classes each of one color. Fermions 
from different classes have mutual semionic statistics. A very intriguing feature of these 
fermions is related to the topological color charges they carry. They carry charges from a 
particular family of low energy fermions. Thus the charges created by open strings are in- 
visible to high energy fermions. Moreover, there are some experimental proposals to realize 
hard-core bosons with optical lattices 91| and it would be a nice challenge to implement a 
Hamiltonian hke (jHj) and ( H5|) . 

We have shown that this new model exhibits enough novel interesting and relevant 
properties so as to justify further research. Some of these possible lines of study are as 
follows: We have only studied a particular phase of the system, although we are able 
to study non-perturbative effects as well. The fact that all phases show a topological 
degeneracy anticipates a rich phase diagram. In this regard, one may explicitly break the 
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color symmetry that the model exhibits and still keep the features that we have discussed. 
It would be particularly interesting to check whether any of the phases displays non-abelian 
anyons. The model has many integrals of motion, although not enough to make it exactly 
solvable. This becomes another appealing feature of the model since other methods of 
study, like numerical simulations and experimental realizations will help to give a complete 
understanding of all its phases. 

Acknowledgements We acknowledge financial support from a PFI grant of EJ-GV, DGS 
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Appendix A: 2-Body Hamiltonian for Color Codes using Cluster States 

A topological color code can be constructed from a graph state defined on a bipartite 
lattice by means of a set of measurements on certain subsystems. This bipartite lattice is 
shown in Fig JTGT a). where the black vertices correspond to plaquettes and the white vertices 
correspond to the vertices of a 2-colex. To this graph we can attach a set of stabilizers as 
follows 

K^ = X^ II Zp, (Al) 

where a and P stand for vertices of the graph and the product runs over all vertices that 
are connected to a by black links. Let us set V = Ui U U2, where Ui and U2 stand for 
the set of white and black vertices of the bipartite graph in Fig JTGT a). Note that white 
and black vertices corresponds to the vertices and plaquettes of the 2-colex. This bipartite 
graph is exactly what we need to construct color codes. To this end, we first impose a 
unitary transformation on the sublattices that allows us to have a more symmetric form of 
the stabilizer operators, i.e 

y feU2 Kf = Z^if), (A2) 

where N{v) denotes the site v and its neighbors, and the same goes for N{f). The corre- 
sponding cluster state denoted by \G) will be the common eigenvector of the above stabilizer 
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(a) 



(b) 



FIG. 16: (color online) (a) The graph needed for obtaining color codes from graph states. The 
graph is bipartite. Black and white vertices correspond to the plaquettes and vertices of a 2-colex. 
Black solid links are edges of the graph (b) The corresponding graph state can be approximated 
as a low energy description of a lattice with 2-body Hamiltonian. The lattice is obtained from the 
graph by replacing its vertices with some hexagons and triangles. The interactions a^a^ and a^a^ 
are associated to the solid and dashed links, respectively. 

operators. Thus we have: 



Finally, a graph state can be related to a color code within a set of measurements in the Z 
basis on all qubits corresponding to the set f/2. 

We suppose there is a two dimensional lattice of physical qubits that is governed by a 
2-body Hamiltonian. Physical qubits of the lattice are projected to logical qubits. The point 
is that this projection is achieved by going to some order in perturbation theory. We think 
of vertices of the graph in Fig JTGT a) as logical qubits. The lattice with 2-body interaction 
is shown in FigIT6l(b). where the number of physical qubits corresponding to the vertices of 
the graph equals the number of links crossing the vertex. The new resulting lattice consists 
of triangles and hexagons and physical qubits live on their vertices. Triangles and hexagons 
are in one to one correspondence with the white and black vertices {Ui and f/2) of the graph 



V / Zn(j)\G) 



\G). 



\G) 



(A3) 
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in Fig JTGT b). respectively. Note that each triangle is linked with three neighboring hexagons 
and each hexagon is linked with six neighboring triangles. 

The low lying spectrum of a well constructed 2-body Hamiltonian defined on the lattice 
composed of hexagons and triangles may describe a cluster state. To this end, we need for 
the following projection from the physical qubits to the logical ones: 

P. = \t.){m\ + i^.)aui, Pf = i^/)(tttttti + \h)miu\ m 

where | ff^) and | with L = v, f stand for the two states of the logical qubits obtained 
within the above projections, or alternatively they are states of logical qubits of the graph 
in FigfT6l(a). We set the following Hamiltonian: 

H = Ho + XV (A5) 

where Hq is the unperturbed Hamiltonian which can be treated exactly and A is a small 
quantity which allow us to treat the term XV perturbatively. We refer to each triangle 
(hexagon) and its vertices by a site index v (/) and indices i, j, respectively. The unperturbed 
part of the Hamiltonian included in Hq is as follows. 

L <i,j> 

where the first sum runs over all triangles and hexagons (sites) and <i,j> stand for the 
nearest-neighbor qubits around the corresponding triangle or hexagon connected by the solid 
lines as in the Fig JIBT b). The interaction between qubits of triangles and qubits of hexagons 
are included in V: 

<vijj> 

where <vi,fj> denotes two neighboring qubits which are connected by the dashed line. 
The unperturbed Hamiltonian Hq can easily be diagonalized for triangles and hexagons. 
Ground state vectors of a triangle or hexagon form a two dimensional space spanned by the 
following polarized eigenvectors: 

triangle: |^,) = |ttt), = (A8) 

and 

hexagon: | ^ | tttttt), \h)^\UUU)- (A9) 
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Let and Nf be the number of triangles and hexagons, respectively. Imposing periodic 
boundary conditions, the total number of triangles and hexagons will be: N = + Nf = 
3Nf . Thus, the dimension of the ground-space of the unperturbed Hamiltonian Hq (or 
ground state degeneracy) becomes: go = 2^"2^f = 2^, and the ground state energy is: 
Eq^^ = Nf{—Q) + A^^(— 3) = —AN in terms of the energy scale of the problem. The first 
excited state is produced by exciting one of the triangles or hexagons and has energy: 
e[^^ = Nf{-6) + A^^,(-3) + 4 = -4(A^ - 1) with degeneracy gi = UN2^ . The second 
excited state has energy E^^^ = Nf{—6) + Ny{—3) + 8 = — 4(A^ — 2) with degeneracy: 
^2 = 4(A^ + 5A^^)2^"\ and so on and so forth. 

Using degenerate perturbation theory as in Sec ilVBl the effect of perturbation V on 
the ground sate subspace can be investigated, and see if it breaks the degeneracy. It is 
simple to see that first order perturbation does not have anything to do with the ground 
state subspace. The second order gives rise to a trivial effect as a shift in energy, since each 
operator related to dashed links appears twice. The third order perturbation theory however 
gives rise to a nontrivial effect. It causes a partial lift of the ground state degeneracy, but not 
complete. The initial degeneracy 2^ gets reduced down to 2'^f . This nontrivial effect arises 
from the product of three dashed links crossing a typical triangle, namely the ground state 
vectors are grouped into the 2^" states, each containing 2^-f vectors. The product of three 
(six) operators around a triangle (hexagon) is equivalent to an X operator acting on the 
logical qubit which is projected down from the three(six) qubits of the triangle (hexagon) , 
since 



Also the action of a on one qubit of a triangle or hexagon is equivalent to an Z operator 
acting on the related logical qubit, since 



XL = \ilL){U\ + \W{ilL 



(AlO) 



ZL = PL(T'Pl=\ilL){ilL 



(All) 



Now we can go on in order to calculate the third order perturbation: 




(A12) 



where 



(A13) 
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the product runs over three black vertices hnked to the v and 5 = §• The operator is a 
stabihzer for the logical qubits which are projected down from the triangles. Since = 1, 
the ground states correspond to the values of k^j = +1. We skip the forth and fifth order of 
perturbation because they have trivial effects. 

Like in S^'^ order perturbation, we are faced with a nontrivial term in the 6^^ order 
perturbation theory. We will see that by considering this order, the ground state degeneracy 
is lifted completely. This nontrivial effect arises from the product of terms in the perturbation 
XV corresponding to the links around a hexagon. Finally, for the Q^^ order perturbation we 
have: 

/ 

where 

Kf = Xfl[Z,, (A15) 

V 

and the product runs over six white vertices linked to the hexagon /. The coefficient 7 has 
positive sign and its precise value is unimportant. We would like to emphasize that at six 
order in perturbation theory some other terms appear which are product of two distinct K^. 
However, we skip them as they all commute. Equations (IA13|) and flA15|) provides all we 
need to adopt the cluster state in (1A3I) as ground state of the low energy effective theory of 



Hamiltonian in (lASp . which up to six order of perturbation can be written as follows 

i^eff = constant — 6 Ky — •y Kf. (A16) 

/ 

We see that the above effective Hamiltonian is completely different from that of in f l34l) . The 
latter equation gives rise directly to the topological color code as its ground state, but the 
ground state (cluster state) of former one needs further local measurements to encode the 
desired color code. 



References 



[1] H. Bombin, M. A. Martin-Delgado; "Topological Quantum Distillation"; Phys. Rev. Lett. 97, 
180501 (2006); Iquant-ph/OeOSTsl 

50 



[2] H. Bombin, M.A. Martin-Delgado; "Topological Computation without Braiding"; Phys. Rev. 

Lett. 98, 160502 (2007); !! quant^h/0610024[ 
[3] H. Bombin, M.A. Martin-Delgado "Exact Topological Quantum Order in D=3 and Beyond: 

Branyons and Brane-Net Condensates"; Phys. Rev. B 75, 075103 (2007); |cond^at/0607736| 
[4] H. J. Briegel, R. Raussendorf, Phys. Rev. Lett. 86, 910 (2001). 

[5] R. Raussendorf, H. J. Briegel, Phys. Rev. Lett. 86, 5188 (2001); Quantum Inf. Comput. 2, 
443 (2002). 

[6] R. Raussendorf, D. E. Browne, and H. J. Briegel, Phys. Rev. A 68, 022312 (2003). 

[7] H. Bombin, M.A. Martin-Delgado, "Statistical mechanical models and topological color 

codes"; Physical Review A 77, 042322 (2008). 'a ?Xrv:0711. 04681 
[8] S. D. Bartlett and T. Rudolph, Phys. Rev. A 74, 040302(R) (2006). 

[9] M. Van den Nest, K. Luttmer, W. Diir and H. J. Briegel, Phys. Rev. A 77, 012301 (2008). 

[10] H. Bombin, M. Kargarian, M.A. Martin-Delgado, "Interacting Anyonic Fermions in a Two- 
Body 'Color Code' Model"; Phys. Rev. B 80, 075111 (20091. 13^0^:0811.0911 1. 

[11] M. Nielsen and I. Chuang, Quantum Computation and Quantum Information (Cambridge 
University Press, Cambridge, England, 2000). 

[12] A. Galindo and M. A. Martin-Delgado, Information and Computation: Classical and Quantum 
Aspects; Rev. Mod. Phys. 74, 347 (2002). Electronic preprint [ arXiv:quant-ph/01 1 21051 

[13] J.E. Levy, A. Ganti, C.A. Phillips, B.R. Hamlet, A.J. Landahl, T.M. Gurrieri, R.D. Carr, 
M.S. Carroll, "The impact of classical electronics constraints on a solid-state logical qubit 
memory" ; |ar2£iY: 0904.0003", 

[14] A. Micheli, G.K. Brennen, P. Zoller, "A toolbox for lattice spin models with polar molecules"; 
Nat. Phys. 2, 341 (2006). |quant-ph/0512222 



[15] L. Jiang, G.K. Brennen, A.V. Gorshkov, K. Hammerer, M. Hafezi, E. Demler, M.D. Lukin, P. 
Zoller, "Anyonic interferometry and protected memories in atomic spin lattices"; Nat. Phys. 



4, 482 f 2008). iarXiv:0711. 13651 
[16] M. Miiller, L. Liang, I. Lesanovsky, P. Zoller, "Trapped Rydberg ions: from spin chains to 

fast quantum gates"; New J. of Phys. 10, 093009 (2008). larXiv: 07 09. 28491 
[17] L.-M. Duan, E. Demler, M.D. Lukin, "Controlling Spin Exchange Interactions of Ultracold 

Atoms in Optical Lattices"; Phys. Rev. Lett. 91, 090402 (2003). arXiv:cond-mat 70210564 



[18] J.J. Garcia-Ripoll, M.A. Martin-Delgado, J.I. Cirac, "Implementation of Spin Hamiltonians 

51 



in Optical Lattices"; Phys. Rev. Lett. 93, 250405 (2004). 
[19] M.J. Hartmann, F.G.S.L. Brandao, M.B. Plenio, "Strongly interacting polaritons in coupled 
arrays of cavities"; Nat. Phys. 2, 849 (2006). |qua^t-ph/0606097l 



[20] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen, U. Sen, "Ultracold atomic gases 
in optical lattices: mimicking condensed matter physics and beyond"; Adv. Phys. 56, 243 
(2007). 

[21] J. Du, J. Zhu, M. Shi, X. Peng, D. Suter, "Experimental observation of a topological phase 
in the maximally entangled state of a pair of qubits"; Phys. Rev. A 76, 042121 (2007). 
larXiv:07 05.3566i 

[22] A.F. Albuquerque, H.G. Katzgraber, M. Troyer, G. Blatter, "Engineering exotic phases for 
topologically protected quantum computation by emulating quantum dimer models"; Phys. 
Rev. B 78, 014503 (2008). 

[23] G.K. Brennen, M. Aguado, J.I. Cirac, "Simulations of quantum double models"; 
larXiv:0901. 13451 

[24] F. Wilczek, "Quantum Mechanics of Fractional-Spin Particles". Phys. Rev. Lett. 49, 957 
(1982). 

[25] J.M. Leinaas, J. Myrheim, "On the theory of identical particles". II Nuovo Cimento B37, 1 
(1977). 

[26] G. Moore, N. Read, "Nonabelions in the fractional quantum Hall effect", Nucl. Phys. B 360, 
362 (1991). 

[27] X.G. Wen, Quantum Field Theory of Many-body Systems: From the Origin of Sound to an 
Origin of Light and Electrons (Oxford Univ. Press, New York, 2004). 

[28] X.G. Wen, "Topological orders and edge excitations in fractional quantum Hall states"; Adv. 
Phys., 44, 405 (1995). 

[29] H. Bombin, M.A. Martin-Delgado, "Interferometry-free protocol for demonstrating topological 

order"; Physical Review B 78, 165128 f2008^. [a?Xrv: 0705.0007[ 
[30] D. Gottesman, "Class of quantum error-correcting codes saturating the quantum Hamming 

bound". Phys. Rev. A 54, 1862 (1996). 
[31] Z.W.E. Evans, A.M. Stephens, "Optimal decoding in fault-tolerant concatenated quantum 

error correction" ; larXiv:0902.4506l 
[32] Z.W.E. Evans, A.M. Stephens, "Accuracy threshold for concatenated error detection in one 

52 



[33 

[34; 
[35; 

[36 

[37; 

[38 
[39 
[40 
[41 

[42; 

[43 
[44; 

[45; 

[46 
[47 



dimension" ; arXiv:0902.2658l 

A. Yu.Kitaev, "Fault-tolerant quantum computation by anyons", Annals of Physics 303, 2 
(2003). |quant-ph/9707021[ 

H. Bombin, M.A. Martin-Delgado; "Homological Error Correction: Classical and Quantum 



Codes"; J. Math. Phys. 48, 052105 (2007); quant-ph/0605094 



H. Bombin, M.A. Martin-Delgado; "Topological quantum error correction with optimal en- 
coding rate"; Physical Review A 73, 062303 (2006). |quant-ph/06020'63| 

CD. Albuquerque, R. Palazzo, E.B. Silva, "Topological quantum codes on compact surfaces 
with genus g>T; J. Math. Phys. 50, 023513 (2009). 

S. Bullock, G.K. Brennen, "Qudit surface codes and gauge theory with finite cyclic groups"; 
J. of Phys. A: Math, and Theor. 40, 3481 (2007). 

H. Bombin, M.A. Martin-Delgado; "Optimal Resources for Topological 2D Stabilizer Codes: 

Comparative Study", Phys. Rev. A 76, 012305 (2007). |quant-p h/0703272 [ 

M. Planat, M.R. Kibler, "Unitary reflection groups for quantum fault tolerance"; 

larXiv:0807.3650l 

L.S. Georgiev, "Computational equivalence of the two inequivalent spinor representations of 
the braid group in the topological quantum computer based on Ising anyons" ; larXiv:0812.2337l 
D. Hu, W. Tang, M. Zhao, Q. Chen, S. Yu, C.H. Oh, "Graphical nonbinary quantum error- 
correcting codes"; Phys. Rev. A 78, 012306 (20081. [a?Xrv:0801.0831[ 

L.S. Georgiev, "Towards a universal set of topologically protected gates for quantum compu- 
tation with Pfaffian qubits"; Nuc. Phys. B 789, 552 (2008). 

H. Bombin, M.A. Martin-Delgado; "Quantum measurements and gates by code deformation" ; 
J. of Phys. A: Math, and Theor. 42, 095302 (2009). larXiv:0704.2540l 

R. Raussendorf, J. Harrington, K. Goyal, "Topological fault-tolerance in cluster state quantum 



computation"; New J. of Phys. 9, 199 (2007). arXiv:quant-ph/0703143 



R. Raussendorf, J. Harrington, "Fault-Tolerant Quantum Computation with High Threshold 



in Two Dimensions"; Phys. Rev. Lett. 98, 190504 (2007). [quant-ph/0610082 

D. P. DiVincenzo, "Fault tolerant architectures for superconducting qubits"; arXiY:Q905- 

E. Dennis, A. Kitaev, A. Landahl, J. Preskill, "Topological quantum memory". J. Math. Phys. 
43, 4452-4505 (2002). 

[48] H.G. Katzgraber, H. Bombin, M.A. Martin-Delgado, "Error Threshold for Color Codes and 

53 



Random 3-Body Ising Models"; Phys. Rev. Lett. 103, 090501 f2009l. [a?Xrv:0902.48451 
[49] M. Ohzeki, "Threshold of topological color code" : larXiv:0903.2102l 

[50] W.-B. Gao, A.G. Fowler, R. Raussendorf, X.-C. Yao, H. Lu, P. Xu, C.-Y. Lu, C.-Z. Peng, Y. 
Deng, Z.-B. Chen, J.-W. Pan; "Experimental demonstration of topological error correction". 
larXiv:0905.1542 V2. 

[51] R. Alicki, M. Fannes, M. Horodecki, "A statistical mechanics view on Kitaev's proposal for 

quantum memories", J. Phys. A: Math. Theor. 40 (2007) 6451-6467. 
[52] R. Alicki, M. Fannes, M. Horodecki, "On thermalization in Kitaev's 2D model"; 

|arXiv:0810.4584l 

[53] R. Alicki, M. Horodecki, P. Horodecki, R. Horodecki, "On thermal stability of topological 

qubit in Kitaev's 4D model" ; larXiv:0811.0033[ 
[54] S. Iblisdir, D. Perez-Garcia, M. Aguado, J. Pachos, "Thermal States of Anyonic Systems"; 

larXiv:0812.4975 i 

[55] S. Iblisdir, D. Perez-Garcia, M. Aguado, J. Pachos, "Scaling law for topologically ordered 

systems at finite temperature" ; larXiv:0806 . 1853 1 
[56] S. Bravyi, B. Terhal, "A no-go theorem for a two-dimensional self-correcting quantum memory 

based on stabihzer codes" ; larXiv:08TO. 19831 
[57] M. Kargarian, "Finite temperature topological order in 2D topological color codes"; Physical 

Review A 80, 012321 f2009). arXiv:0904.4492 . 
[58] M. Kargarian, "Entanglement properties of topological color codes"; Physical Review A 78, 

062312 (2008). larXiv:0809.4276l 
[59] H-H. Tu, G-M. Zhang, T. Xiang, Z-X. Liu, T-K. Ng, "Topologically distinct classes of valence 

bond solid states with their parent Hamiltonians" ; larXiv : 0904jO_55D 1 
[60] A. Pakman, A. Parnachev, "Topological entanglement entropy and holography"; JHEP 07, 

097 (2008). arXiv:0805.189n 
[61] E. Rico, H.J. Briegel, "2D multipartite valence bond states in quantum anti-ferromagnets" ; 

Ann. of Phys. 323, 2115 (2008). arXiv:0710.2349l 
[62] J. Eisert, M. Cramer, M.B. Plenio, "Area laws for the entanglement entropy - a review"; 

l arXiv:0808.3773[ 

[63] W. Diir, H.J. Briegel, "Entanglement purification and quantum error correction"; Rep. on 
Prog, in Phys. 70, 1381 (2007). arXiv:0 705.4165l 

54 



[64] V. Karimipour, "A complete characterization of the spectrum of the Kitaev model on spin 



ladders" ; larXiv:0904.3554l 

[65] M. Van den Nest, W. Diir, and H. J. Briegel, "Classical Spin Models and the Quantum- 
Stabilizer Formalism"; Phys. Rev. Lett. 98, 117207 (2007); | arXiv:quant-ph/0610157] 

[66] S. Bravyi and R. Raussendorf, "Measurement-based quantum computation with the toric code 



states"; Phys. Rev. A 76, 022304 (2007); |arXiv:quant-ph/0610162f 



[67] M. Van den Nest, W. Diir, and H. J. Briegel, "Completeness of the classical 2D Ising model 

and universal quantum computation" ; larXiv:0708.2275l 
[68] M. Van den Nest, W. Diir, R. Raussendorf, H.J. Briegel, "Quantum algorithms for spin models 

and simulable gate sets for quantum computation" ; larXiv:0805.1214l 
[69] R. Hiibener, M. Van den Nest, W. Diir, H.J. Briegel, "Classical spin systems and the quantum 

stabilizer formalism: general mappings and applications" ; larXiv:08i2.2l27l 
[70] G. De las Cuevas, W. Diir, M. Van den Nest, H.J. Briegel, "Completeness of classical spin 

models and universal quantum computation"; larXiv:0812.23 68 . 
[71] G. De las Cuevas, W. Diir, H. J. Briegel, M. A. Martin-Delgado, "Unifying all classical spin 

models in a Lattice Gauge Theory"; Phys. Rev. Lett. 102, 230502 f2009). iarXiv:0812.3583l 
[72] A. Yu. Kitaev; "Anyons in an exactly solved model and beyond"; Ann. of Phys. 321, 2 (2006). 

|arXiv:cond-mat /0506438V3. 



[73] X.-Y. Feng, G.-M. Zhang, T. Xiang, "Topological Characterization of Quantum Phase Tran- 



sitions in a Spin-1/2 Model"; Phys. Rev. Lett. 98, 087204 (2007). [coml-mat/0610626 
[74] H. Yao, S.A. Kivelson; "Exact Chiral Spin Liquid with Non-Abelian Anyons"; 

Phys. Rev. Lett. 99, 247203 (2007). larXiv:0708.0040l 
[75] G. Baskaran, S. Mandal, R. Shankar, "Exact Results for Spin Dynamics and Fractionalization 

in the Kitaev Model"; Phys. Rev. Lett. 98, 247201 (2007). 
[76] S. Yang, D.L. Zhou, CP. Sun; "Mosaic spin models with topological order"; Phys. Rev. B76, 



180404 (R) (2007). |ax2Oy:0708.0676l 
[77] T. Si, Y. Yu, "Anyonic loops in three-dimensional spin liquid and chiral spin liquid"; Nuc. 

Phys. B 803, 428 (2008). [a?Xiv:0712.4231[ 
[78] S. Mandal, N. Surendran, "Exactly solvable Kitaev model in three dimensions"; Phys. Rev. 

B 79, 024426 (2009). 

[79] M. Levin, X.-G. Wen, "String-net condensation: A physical mechanism for topological phases 

55 



"; Physical Review B71, 045110 (2005). 
[80] M. Oshikawa and T. Senthil, "Fractionalization, topological order, and quasiparticle statistics"; 

Phys. Rev. Lett. 96, 060601 (2006). 
[81] D. L. Bergman, R. Shindou, G. A. Fiete, and L. Balents, "Degenerate perturbation theory of 

quantum fluctuations in a pyrochlore antiferromagnet" ;Physical Review B 75, 094403 (2007). 
[82] G. Kells, A. T. Bolukbasi, V. Lahtinen, J. K. Slingerland, J. K. Pachos, and J. 

Vala, "Topological Degeneracy and Vortex Manipulation in KitaevOs Honeycomb Model"; 

Phys. Rev. Lett. 101, 240404 (2008). 
[83] C. Knetter, G. Uhrig; "Perturbation theory by flow equations: dimerized and frustrated S = 

1/2 chain"; Eur. Phys. J. B13, 209 (2000). 
[84] F. Wegner, "Flow-equations for Hamiltonians" ; Ann. der Phys. 506, 77 (1994). 
[85] J. Vidal, K.P. Schmidt, S. Dusuel; "Perturbative approach to an exactly solved problem: 

Kitaev honeycomb model"; Phys. Rev. B 78, 245121 (2008). arXiv:0809. 1 5531 
[86] M. Levin, X.-G. Wen, "Fermions, strings, and gauge fields in lattice spin models"; Physical 

Review B 67, 245316 (2003). 
[87] H. Bombin, M.A. Martin-Delgado, "Family of non-Abelian Kitaev models on a lattice: Topo- 
logical condensation and confinement"; Physical Review B 78, 115421 (2008). 
[88] H. Bombin, M.A. Martin-Delgado, "Nested Topological Order" : [arXiv:0803.4299 l 
[89] F.A Bais, J.K. Slingerland, "Condensate-induced transitions between topologically ordered 

phases"; Physical Review B 79, 045316 f2009). [aFxrv:0808.0627l 
[90] Ehiott H. Lieb, "The flux phase of the half-filled band"; Phys. Rev. Lett. 73, 2158 (2009). 
[91] B. Capogrosso-Sansone, C. Trefzger, M. Lewenstein, P. Zoller, G. Pupillo; "Quantum Phases 

of Cold Polar Molecules in 2D Optical Lattices" : [arXiv:0906.2009l 



56 



